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Abstract. 

We investigate the field emission properties of nitrogenated and boronated carbon 
nanotubes using time-dependent density functional theory, were the wave function 
propagation is performed using the Crank-Nicholson algorithm. We extract the 
current-voltage characteristics of the emitted electrons from nanotubes with different 
doping configurations. We found that boron doping either impedes, or slightly 
enhances, field emission. Nitrogen strongly influences the emission current, and the 
current is sensitive to the location of the nitrogen dopant in the nanotubc. The emitted 
charge cloud from nitrogen doped carbon nanotubes is, however, more diffuse than that 
from pristine ones, our simulations show the emergence of a branching from the charge 
cloud, making nitrogenated carbon nanotubes less convenient for use in narrow beam 
applications. 
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1. Introduction 

Field emission has been one of the active research areas owing to its theoretical, as 
well as commercial, significance pLj. Unique field emission properties were discovered in 
carbon nanotubes owing to their high aspect ratio and mechanical and chemical stability, 
making them the best candidate for flat-panel displays and atomic force microscopy 
probes [21 El S] - Among the number of nanostructures investigated for field emission, 
carbon nanotubes proved, both theoretically and experimentally, to possess strongest 
field emission characteristics. 

Nitrogen and boron are among the most common dopants encountered in carbon 
nanotubes [5]. Nitrogen-doped carbon nanotubes were studied experimentally by several 
groups [6]|7][8][9], reporting generally improved field emission characteristics. However, 
the strength of the emission current, as well as the onset field, were reported to 
be influenced by the substrate used in fabricating the nitrogenated carbon nanotube 
[S], processing method [8] and post-processing [7]. While nitrogen doping of carbon 
nanotubes, as was indicated by, for example, Doytcheva et al. [8], is close to that 
of metallic emitters, boron doping either impedes emission or has minimal effect on 
emission. Sharma et al. \IU\ suggested that, based on their experimental investigation, 
boron doped carbon nanotubes are more convenient for applications in flat panel display 
devices where they are packed in large arrays or parallel emitters, while nitrogen-doped 
carbon nanotubes, owing to their high field emission quality, can be used in narrow 
beam applications. 

Theoretical work by Qiao et al. [11], and experimental work by Chan et al. [12] 
concluded that boron doping decreases the tunneling probability in comparison with 
pristine carbon nanotubes. Owens [15] . based on theoretical computations, expected 
that both boron and nitrogen doping will not enhance field emission based on the fact 
that they increase the ionization energy of carbon nanotubes. However, although this 
is true in the case of boron doping, the ionization potential was not a correct indicator 
for field emission quality. Theoretical work on the field emission from boron doped 
carbon nanotubes, however, have reported conflicting results |H] , and there is a need to 
understand the details of the emission process and the variable factors to which emission 
current in boron doped carbon nanotubes would be sensitive. We believe that this work 
will shed more light into the microscopic nature of field emission from doped carbon 
nanotubes and the role of electron donation or acceptance in enhancing or impeding 
field emission. 

2. Computational Details 

We summarize our calculation procedure as follows (more details of the procedure are 
presented in Ref. ??). We performed Density Functional Theory computation using 
the Local Density Approximation for exchange and correlation potential by Pedrew and 
Zunger [14J. Then, we solved the time-dependent Kohn-Sham equations by propagating 
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the wave function using the Crank-Nicholson method. Then, following each propagation 
step, the amount of charge remaining in the emitter region at time t for a particular 
wave function corresponds to the electric charge remaining, which is given by 

0»(*) = jf7 / H 2rfx > C 1 ) 
and in terms of the life time r n of state n, 

Q n (t) = e-%, (2) 

which implies linear behavior of Q n (t) in the short time interval. From this formula, 
the current generated by a wave function is given by 

/. = e ^„-l (3) 

in the short time range. The total current is given by 

^ J dt 

In order to compute the charge remaining in the lower portion of the simulation box, 
we integrate that lower portion each time-step for each wave function. At the beginning 
of the simulation, the charge is typically 1 (which works as a check for normalization 
consistency of the propagator used). Charge starts to decrease as time elapses, until it 
reaches a minimum point after which it starts increasing (due to reflection of the wave 
back from the upper surface of the simulation box). The integration is performed by 
simply summing the product of unit volume boxes of the mesh used in Octpus within 
the region concerned. 



3. Results and Discussion 



Here we consider three different nitrogen and boron doping configurations: the dopant 
replacing a carbon atom in the top ring; the dopant replacing a carbon atom in the 
second ring; the dopant replacing a carbon atom in the third ring. The different positions 
are shown in Fig. [T] below. 

In order to understand the mechanism whereby a dopant enhances or impedes 
emission, we present in Fig. ?? the charge evolution graphs for both the Nl and B2 
configurations. 

Fig. [2] indicates the occurrence of a strong propagation of the wave function into 
the vacuum region. The influence of the nitrogen atom, however, did not concentrate 
the charge into a specific direction. Charge dispersed into a region that is almost 
wide. Comparison with the charge evolution graph for pristine carbon in Fig. [4] shows 
the difference between the structure of the emitted charge cloud in pristine and nitrogen 
doped carbon nanotubes. In the former, the charge cloud tends to concentrate at specific 
regions, or channels, whereas in the latter case, the cloud tends to disperse in space. 
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Figure 1. Relative positions of nitrogen (or boron) dopant on the carbon nanotube 
tip 

Table 1. Extracted current (/xA) from prestine vs. nitrogen-doped/boron-doped 
carbon nanotube against an applied electric field of 0.2V/ A. All are H-passivated with 
no extra electrons. Nl: Nitrogen dopant replacing a carbon atom at the top ring. 
N2: Nitrogen dopant replacing a carbon atom at the second ring below the tip. N3: 
Nitrogen dopant replacing a carbon atom at the third ring below the tip. Bl: Boron 
dopant replacing a carbon atom at the top ring. B2: Boron dopant replacing a carbon 
atom at the second ring below the tip. B3: Boron dopant replacing a carbon atom at 
the third ring below the tip 



Configuration 


Current 


Pure 


0.0161 


Nl 


0.53394 


N2 


0.0796 


N3 


0.0164 


Bl 


0.02198 


B2 


0.013469 


B3 


0.0268016 



It is an interesting observation that in nitrogen doping, the emission current is 
highly sensitive to the location of the dopant, which is not the case in boron doping. 
The location of a nitrogen atom at the uppermost ring (the tube's apex) results in a 
greatly enhanced emission current. As we place the nitrogen dopant further away from 
the tip, the emission enhancement drops drastically. This observation can be explained 
on the basis of the presence of the donated electron away from the tube body and closer 
to the tip region which already has increased electron localization. At the tube's apex, 
the donated electron enriches the heavy charge cloud in the tip region, and being away 
from the body of the tube, it experiences lower ionization potential. As the donated 
electron is placed further away from the apex, it experiences higher ionization potential 
imposed by the tube's body. 
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Figure 2. Evolution of the wave function into the vacuum region in nitrogen-doped 
carbon nanotube: highest energy wave function. Units are arbitrary. 

As for boron doping, results in Tab. [T] indicate that a boron dopant placed at the 
second top ring in the nanotube will impede emission, whereas placing it either at the 
apex or further down will enhance emission by increasing the emission current by 25% 
compared to emission from a pristine carbon nanotube. Fig. [3] below shows the impact 
of the boron dopant: the charge cloud seems to be almost fixed. This is due to the fact 
that boron is a charge acceptor; it attracts the electron cloud, contrary to the case of 
nitrogen doping. 
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Figure 3. Evolution of the wave function into the vacuum region in boron-doped 
carbon nanotube: highest energy wave function. Units are arbitrary. 




Figure 4. Evolution of the wave function into the vacuum region in pristine carbon 
nanotube: highest energy wave function. Units are arbitrary. 

Fig. [5] shows the difference between the charge cloud emerging in pristine, boron- 
doped and nitrogen-doped carbon nanotubes subject to 0.8V/A, where the dopant is 
placed at the nanotube's apex. Each column shows the evolution of field emission from 



First Principles Calculation of Field Emission 



7 



the tip region of the nanotube after 10 au, 20 au and 30 au. The boron-doped carbon 
nanotube hardly emits charge. The extent of the charge cloud along the x axis in boron 
is smaller than that in pristine, which is yet smaller than that in nitrogen-doped CNT. 
This is a consequence of the fact that boron is a charge acceptor, whereas nitrogen is 
a charge donor. Nitrogen doping clearly enhances the field emission by expanding the 
size of the charge cloud surrounding the tip region, and by liberating more charge as 
time elapses (which is depicted in the figure as a brightening of the blue color). 




Figure 5. Evolution of the wave function into the vacuum region in nitrogen-doped 
carbon nanotube: highest energy wave function. Units are arbitrary. 

The extent to which the emerging charge cloud is diffuse in nitrogen-doped carbon 
nanotube can be inferred from Fig. [6] below. Whereas the pristine carbon nanotube 
shows the emergence of a well-defined protrusion directed upwards, the charge cloud 
emerging from the nitrogen-doped nanotube diffuses to the surrounding tubes. Thus, 
although the signal from the nitrogen-doped nanotube is much stronger, it will be much 
harder to attempt to project it towards a particular target. 
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Figure 6. Evolution of the wave function into the vacuum region in nitrogen-doped 
carbon nanotube: highest energy wave function. Units are arbitrary. 

4. Conclusion 

In summary, we have performed time-dependent density functional theory simulation 
of the field emission from nitrogen- and boron- doped carbon nanotubes. Boron tends 
to either weaken or slightly improve the emission current whereas nitrogen strengthens 
it. The strength of the emitted current in nitrogen-doped carbon nanotube is highly 
sensitive to the position of the dopant. And although the current emerging from the 
nitrogen-doped carbon nanotube is much stronger than that of pristine carbon nanotube, 
the charge cloud emerging from the latter is seen to project into a specific direction, 
unlike that of the former, which is diffused into space. We believe that the consideration 
raised in this paper about the extent of the sharpness of the emerging charge cloud will 
be of utility to applications of carbon nanotubes involving narrow electron beams, such 
as in Tunneling Electron Microscopy. Much work still needs to be done in identifying 
the dopants, as well as the conditions, that can produce narrower charge clouds. 
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